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Adiabatic quantum annealing is a paradigm of analog quantum computation, 
where a given computational job is converted to the task of finding the global mini- 
mum of some classical potential energy function and the search for the global poten- 
tial minimum is performed by employing external kinetic quantum fluctuations and 
subsequent slow reduction (annealing) of them. In this method, the entire potential 
energy landscape (PEL) may be accessed simultaneously through a delocalized wave- 
function, in contrast to a classical search, where the searcher has to visit different 
points in the landscape (i.e., individual classical configurations) sequentially. Thus 
in such searches, the role of the potential energy might be significantly different in 
the two cases. Here we discuss this in the context of searching of a single isolated 
hole (potential minimum) in a golf-course type gradient free PEL. We show, that 
the quantum particle would be able to locate the hole faster if the hole is deeper, 
while the classical particle of course would have no scope to exploit the depth of the 
hole. We also discuss the effect of the underlying quantum phase transition on the 
adiabatic dynamics. 



I. INTRODUCTION 

Adiabatic quantum annealing (AQA) is a method of finding the ground state 

(minimum energy state) of a given classical Hamiltonian by employing external quantum 
fluctuations and subsequent adiabatic reduction of them. One is given with a classical 
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Hamiltonian Ti., which may be a physical Hamiltonian with many degrees of freedom, or a 
suitable mathematical function depending on many variables, and the task is to determine 
its global minimum. In order to introduce the quantum fluctuations necessary for the AQA 
of such a Hamiltonian, one adds a quantum kinetic part TC'(t) to it, such that TC'(t) and 
7i do not commute. Initially, one keeps \H'(t — 0)| ^> \H\ so that the total Hamiltonian 
7~Ltot(t) = T~t'(t) + H is well approximated by the kinetic part only (7Y iot (0) ~ 7Y'(0)). If 
the system is initially prepared to be in the ground state of H'(0) (one chooses H'(0) to 
have a easily realizable ground state) and TC'it) is reduced slowly enough, then according 
to the adiabatic theorem of quantum mechanics, the overlap \(ip{t)\E^(t))\, between the 
instantaneous lowest-eigenvalue state \E_(t)) and the instantaneous state of \ip(t)) of the 
evolving system, will always stay near its initial value (which is close to unity, since 'Htot{^>) ~ 
Tl'(p)). Hence at the end of such an evolution, when H,'(t) is reduced to zero at t — r (the 
annealing time), the system will be found in a state \iP{t)) with \(ip(r)\E_(r))\ ~ 1, where 
|£L(t)) is the ground state of Titotij), which is nothing but the surviving classical part TC. 
Thus at the end of an adiabatic annealing the system is found in the ground state of the 
classical Hamiltonian with a high probability. Based on this principle, algorithms can be 
framed to anneal complex physical systems like spin glasses as well as the objective functions 
of hard combinatorial optimization problems (like the Traveling Salesman Problem), towards 
their ground (optimal) states 
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In order to ensure adiabaticity, the evolution should be such that 
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\E + (t)) being the instantaneous first excited state of the total Hamiltonian 7itot{t), A(t) 
being the instantaneous gap between the ground state and the first excited state energies 
and a being the adiabatic factor (for a simple proof see 11410. 
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One key feature, believed to be behind the success of AQA over the classical ones [2J, [6|, [8( 
in glass-like rugged PEL, is the ability of the quantum systems to tunnel easily through 
potential energy barriers even if they are very high, provided they are narrow enough, in 



3 



contrast to the the classical ones, which always has to scale the barrier height with its kinetic 
energy (temperature) irrespective of the width {|], ?], Q, Q, Q- Here we show that there 
is another aspect which makes a quantum mechanical searcher more advantagous over the 
classical ones- it can utilize the depth of the potential energy minimum in locating it in 
absence of any potential gradient which a classical searcher cannot. 

II. SEARCHING A HOLE ON A GRADIENT-FREE PEL 

We consider a lattice with N sites, \i) denoting the state of a particle localized at the 
z-th site. At each site, there is a potential, which is zero at all the sites i ^ w, and is — \ 
at % — w, where w is chosen randomly. Thus the PEL is essentially a flat one without any 
gradient, with a single hole (minimum) at i = w with a depth \. This is precisely some kind 
of analog version of Grover's algorithm for searching a particular entry in an unstructured 
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database [16j - [18[ . But in those studies, the possibility of utilizing the depth of the hole in 
favor of faster search was not considered, and the gain over the classical ones is limited by 



the optimal Grover's bound of 0{y ~N) speed-up [18l. Il9l|. 

Let us consider that the lattice points are connected to each other by an infinite range 
hopping term T between any two sites. The question is how fast a particle can locate the hole 
starting from a state which does not assume any knowledge of the position of the hole, by 
reducing its kinetic energy T, and tuning the hole depth x- The Hamiltonian for a particle 
on such a lattice will be given by- 

nut) = - x (t)\w)(w\-r(t) KM; x{t),m > 0. (3) 

In order to anneal the particle to the hole, one has to reduce T from a very high value to 
a very low final value and tune x m the opposite manner, so that F(t = 0) 3> x(t = 0) 
and T(t = r) x{t = T )i where r is the annealing time. The evolution should satisfy the 
adiabatic condition (TjQ). The eigen-spectrum of 7i tot (t) consists of a ground state \E_(t)) 
and first excited state \E + (t)) (in the order of increasing eigen values) with energies 

E±(t) = - 1 - [(JV - 2)r + x ± V(NT - x) 2 + 4*r] (4) 

respectively, all the time dependencies being implicit, through the time dependence of T and 
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X- The instantaneous gap is thus given by 



A(t) = v/(AT- X )2 + 4 X r . 



(5) 



The instantaneous first and second excited states \E±(t)) are given by 
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where 



C ± (t) = ^[-(N-2)T + XT A] 



(7) 



The second excited state is (N — 2)-fold degenerate, with eigenvalue — r, and the time 
evolving Hamiltonian never mixes the first two eigenstate with any of the second excited 
states. This can be easily argued noting that a state of the form j^) = ^(K) — \j)) 
7^ w) is an eigenstate of H.tot{t) with eigenvalue — T, and (E2\E-(t)) = (E2\E + (i)) = 
for all t. For all allowed combinations of i and j we get (N — 2) such linearly-independent 
eigenstates. Form these (N — 2) linearly-independent eigenstates we can construct (N — 2) 
mutually orthogonal eigenstates, each of which will obviously satisfy the above non-mixing 
condition. Thus we have to take care of only two lowest lying states and the gap between 
them. Henceforth we will consider only the large N limits, and replace ' =' by ' rs' when 
the correction will vanish in the said limit. 

In investigating the condition of adiabatic evolution we consider two separate cases. In 
general, in a quantum annealing program, one might not have the facility of tuning both 
the potential energy part and the kinetic part in practice, since, say, one might depend 
on the strength of interactions between the elementary constituents (hard to tune), while 
the other might be introduced through an applied external field (easily tunable). Hence in 
our analysis we have considered two separate cases - in the first case we tune x keeping V 
constant, while in the second case we do the reverse. 



We start with the initial state |-0(O)) = ^= £V \i) (which is of course the ground state 
of T~ttot(ty) an d adopt a linear annealing schedule: 



A. Constant-r Annealing 



r = r o; x = xo-; xo = rivr . 
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FIG. 1: Panel (a) in the figure shows the variation of the minimum time T m in required to achieve the target 
success probability Pt — 0.33 for the constant-r annealing (T(t) = T , x(t) = X^I T \ r — 2.0, r = 0.5), 
obtained by solving the time-dependent Schrodingcr equation numerically. Panel (b) shows the variation of 
final probability P(\w}) of finding the system in state \w), with annealing time r, for different final value of 
r. Here we have taken N = 10 6 . 

Here r is an arbitrary factor which might also depend on N. When T is kept constant, 
we have to keep it sufficiently low (or, in other words, xo sufficiently high), so that the 
ground state of TCtot(T~) has a substantial (non-vanishing in the N — > oo limit) overlap with 
\w). At t = r, we get from Eq. ([7]), the amplitude of \w) in the final ground state to be 
C-{t)/^C 2 (t) + N - 1, where 
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Clearly, if r > 1, then the amplitude C_(r) > N, and thus the overlap amplitude ~ 0(1), 
whereas if r < 1, then C_(r) ~ 1, and the amplitude vanishes as iV — > oo. Thus to be able 
to locate the hole at the end, we have to take r > 1. In fact, The adiabatic factor in that 
case (r > 1) is given by 

\{nUt))\ 
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This has its maximum at t m = r(N — 2)r /xo ~ T l T i an d the maximum value is given by 
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This clearly shows that if the hole-depth xo scales linearly with the size N of the search 
space (i.e., r is independent of N), the search can be completed in a time independent of N. 
We calculated numerically the minimum time r m j„ required for obtaining a targeted succes 
probability Pt = \ (iP(t)\w)\ 2 = 0.33 for different N, through many decades. The evolution 
is computed solving time-dependent Schrodinger equation numerically and r min is figured 
out up to an accuracy of 10~ 4 employing the following bisection scheme. We first figure out 
arbitrarily a high value of r (call it 77^) for which P(\w)) = \(i/j(t)\w)\ 2 > Pt- Next we find 
a low r (rio), for which P(\w)) < Pt- Then we evaluate P(\w)) for r = r m = {r hi + ti q )/2. If 
the result is greater than Pt, then we replace rjy by r m (and retain old ti ), else we replace 
rio by r m (and retain old Tjy) and repeat the same process. We go on iterating until the the 
value of \P(\w))—Pt\ for both tm and r/ Q lies within some desired accuracy limit. The results 
(Fig. [Tk) clearly show that P(\w)) becomes independent of N for large N, as expected. 

The relaxation behavior for large N for a given annealing time r of course depends on 
the value of T (see Fig. [TJd). If r is too small, the system takes a longer time to feel 
the changes in the landscape, and hence the adiabatic relaxation requires longer time (the 
adiabatic factor becomes bigger; see Eq. [TTil . On the other hand, if To is too large, the 
ground state itself is pretty delocalized, and hence the final state, though more closer to the 
ground state, has again a small overlap with the target state \w). For T « 0.5, the schedule 
is found to be optimal (Fig. QJb)). The relaxation behavior is seen to be linear with the 
annealing time r for large r. 

B. Constant-x Annealing 

Next, we consider the case where x is kept fixed and T is reduced linearly, i.e., x = 
Xo = rNTo and T = Tq(1 — t/r). In this case, we start with the same democratic initial 
state j^o) as in the previous case, which is not of course the ground state in presence of the 
hole (we cannot construct the actual initial ground state without the explicit knowledge of 
the location of the hole). All we need to show in this case, is that it has a non- vanishing 
(i.e., non-zero in the N —>■ 00 limit) overlap |(^(0)|£'_(0))| between our initial state and 
the true ground state at t — 0. If we can assure adiabaticity for the subsequent evolution, 
this overlap will be conserved and would emerge as the final overlap between \ip{r)) and 
\w). Calculation similar to that of the previous section shows that, for r < 1, C_(0) — > as 
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FIG. 2: The figure shows variation of the minimum time T min required to achieve success probability 
Pt = 0.9, with TV for the constant-x annealing. The result is obtained by solving the time-dependent 
Schrodinger equation numerically with T(t) = Tq(1 — t/r), x(t) = Xo = rToN for r = 0.5. 

N — > oo, implying that the overlap actually tends to unity in that limit, while it vanishes 
for r > 1. The adiabatic factor in the former case (r < 1) is given by 

1 



a 2 
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which is again a constant independent of N, if r is iV-independent (i.e., Xo is linear in N). 
We present a similar numerical calculation, as in the previous section, for the variation of 
Tmin with N in Fig. ([2]). It shows that the value of r m « n becomes independent of N in the 
large- N limit (its asymptotic value, however depends on r ). 



III. THE UNDERLYING TRANSITION 



Let us concentrate on the rather interesting case of iV-independent r. The expression 
for the gap A (Eq. [5]) shows that the gap scales linearly with N everywhere except at the 
point iVF = T c = x, where it scales as yN. In our notation, this happens at the instant 
t* — t/t in the case of constant-r annealing, and at t** = r(l — r) in the constant-x case. 
For successful annealing one cannot avoid this point during the evolution in both cases. To 
see what happens at that point, let us focus on the amplitude \(w\E_(x,T)) for T > r c 
and T < r c . Let us consider, say, the T-constant annealing case. We have from Eq. (j7j), 
C-(t) ~ (— iVTo(l — rtjr) + iVTo|l — rt/r\) /2Yq. Taking the modulus into account, one finds 
C-(t) ~ for t < t* (which means the amplitude vanishes), while C_(t) ~ for t > t* (the 
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amplitude tends to unity). This means the ground state of H-totix^) undergoes a global 
change in character from a completely delocalized one to a completely localized one at this 
special point. Thus, to follow the ground state, a lot of tunneling from all the sites to the 
hole is required at the point. Here the depth of the hole plays a crucial role in making this 
tunneling possible preventing the gap A from closing at the transition point. If however, 
the system passes this point very fast, this massive tunneling might remain incomplete and 
a resultant loss in the adiabaticity may occur. Similar argument holds for the x- cons t an t 
annealing. 

If instead of the Hamiltonian we considered, one takes a bounded version of it (energy 
not growing with N), say, by keeping \ independent of N and scaling Tq by N, then one 
can easily see that this special transition point would behave like a true quantum critical 
point with the gap vanishing as 1/yN 17]. The two states- the one localized at the hole 
and the other delocalized over the entire lattice, have exactly equal energies at this "critical" 
point as the gap closes. But before reaching that point (i.e., in the kinetic energy dominated 
region), the ground state is non-degenerate and delocalized, while after the point (i.e., in 
the potential energy dominated region) the ground state is degenerate and localized at the 
hole. Since there is no energy difference between the two states (the localized and the 
delocalized ones) at the critical point, the evolving system cannot sense the global exchange 
of character between the ground state and the first excited state at this point and thus fails 
to adopt (tunnel) accordingly to follow the ground state. Letting the hole scale linearly with 
iV and taking infinite range hopping, one prevents the gap from closing, and thus favors the 
localized state energetically over the delocalized one at the transition point. This acts as a 
drive for the necessary tunneling and thus provides a guidance for the system to follow the 
ground state as it changes its global character. 

However, it is also worth noting that increasing the hole-depth xo indiscriminately does 
not pay. In the case of constant-r annealing, the adiabatic factor increases with r (Eq JTTj) . 
which means we have to keep r as small as possible, i.e., r — > 1 gives the best result. This 
can be explained by noting that a larger xo would demand a slower rate of evolution if 
adiabaticity is to be ensured. In the constant-x case, adiabatic factor decreases with r, but 
one has to keep r < 1, which means x can be increased only linearly with N, up to the 
upper limit NT in order to accelerate the annealing effectively. This is because, too big Xo 
implies too much error in the initial state and the overlap is negligible right from the onset. 
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Hence even with a perfectly adiabatic evolution we end up with the negligible overlap. 

To summarize, we studied the adiabatic searching of a gradient-free potential energy 
landscape on a lattice with a randomly placed single isolated minimum (hole) by a quantum 
searcher. We have found that the depth of the hole plays a non-classical role in accelerating 
the search. This is because the delocalized wave-function of the quantum adiabatic searcher 
can detect the hole from the very onset of the search and keep track of it, in contrast to a 
classical searcher, who of course cannot sense the hole depth non-locally and utilize it. We 
found that the condition for adiabatic search requires an infinite-range hopping between the 
points in the PEL and the depth of the hole to scale linearly with the lattice-size N. We 
have discussed how the hole-depth plays a role in preventing the closure of the gap and the 
consequent divergence of the searching time. 
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